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ABSTRACT 

Binary mergers (NSMs) of double neutron star (and black hole-neutron star) systems are suggested 
to be major sites of r-process elements in the Galaxy by recent hydrodynamical and nucleosynthesis 
studies. It has been pointed out, however, that the estimated long lifetimes of neutron star binaries are 
in conflict with the presence of r-process-enhanced halo stars at metallicities as low as [Fe/H] ~ —3. 

To resolve this problem, we examine the role of NSMs in the early Galactic chemical evolution on the 
assumption that the Galactic halo was formed from merging sub-halos. We present simple models for 
the chemical evolution of sub-halos with total final stellar masses between Mq and 2 x 10® Mq. 
Typical lifetimes of compact binaries are assumed to be 100 Myr (for 95% of their population) and 
1 Myr (for 5%), according to recent binary population synthesis studies. The resulting metallcities 
of sub-halos and their ensemble are consistent with the observed mass-metallicity relation of dwarf 
galaxies in the Local Group, and the metallicity distribution of the Galactic halo, respectively. We 
find that the r-process abundance ratios [r/Fe] start increasing at [Fe/H] < —3 if the star formation 
efficiencies are smaller for less massive sub-halos. In addition, the sub-solar [r/Fe] values (observed as 
[Ba/Fe] ^ —1.5 for [Fe/H] < —3) are explained by the contribution from the short-lived (~ 1 Myr) 
binaries. Our results indicate that NSMs may have a substantial contribution to the r-process element 
abundances throughout the Galactic history. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — Galaxy: evolution — Galaxy: 
halo — galaxies: dwarf — stars: abundances — stars: neutron 


1. INTRODUCTION 

The astrophysical origin of the rapid neutron-capture 
(r-process) elements is still unknown. Recent spectro¬ 
scopic studies of Galactic halo stars have revealed the 
presence of Eu (almost pure r-process element in the 
solar system) in the atmospheres of extremely metal- 
poor (EMP) stars down to [Fe/H] 0 —3 with a 

large star-to-star scat ter in [Eu/Fe] (more than 2 or¬ 
ders of magnitude, e . g..lHonda et al.l20M : iFrancois et al.l 
l2007t iSneden et al.l 1200^ These observational fea¬ 
tures are reasonably reproduced by models of inhomoge¬ 
neous Galactic chemical evolution (GCE), where a lim¬ 
ited mass range (e.g., 8 — IOMq) of core-collapse super¬ 
novae (CCSNe) is assumed to be the dominant source 
of the r-process elements lllshimaru &: Wanaid 119991 

Tsuiimo to. Shige vama. fc YoshiirT200rt llshimaru et al l 

2004 lArgast et al.ll2004H . Previous homogeneous GGE 
models that assume a well-mixed, one-zone Galactic halo 
also suggest the low-mass end of CGSNe as the major r- 
process site ([Mathews et al.l[T992HTravaglio et al.lli999ll 
but are not capable of accounting for the star-to-star 
scatter in [Eu/Fe]. 

To date, however, no study of CCSN nucleosynthesis 
satisfactorily accounts for the production of heavy r- 
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[A/B] = \o^{Nj\/Nb) — \q^{Nj\,^/Nb)q‘> where Nj\ indicates 
the abundance of A. 


proce ss elements (see, e.g.. iWanaio et aOl20lH IWanaiol 
IMI). As an alternative scenario, binary mergers of dou¬ 
ble neutron stars (NSMs) as well as of neutron star-black 


hole pairs have long been su 

ggested to be the r-process 

sites ( 

Lattimer & SchrammI 

1974 119761: ILattimer et al.l 


Svmbalistv & SchrammI 119821: Eichler et al.l 119891: 

Mever 

studie 
NSMs 
in the 

1198911. In fact, a number of recent nucleosynthesis 
s based on hydrodynamical simulations support 
to be the promising sources of r-process elements 
Galaxv (iFreiburghaus et al.l 119991: IGorielv et al.l 

20111: iKorobkin et al.l 120121 Bauswein et al.l 120131: 


Rosswog et al.ll20l4 IWanaio et al.ll2014ll . 


Recent GCE models appear to disfavor the NSM 
scenario except for the cases invoking ve ry short life- 


times of neutron star binaries fl-10 Myr. lArgast et al. 

2004 

: IDe Ponder & VanbeverenI 12004 IMatteucci et al. 

2014 

: IKomiva et al. 

[2014 

Tsuiimoto & Shigevamal 

2014 

van de Voort et al.l 

2015|: 

Wehmever et al.l 1201511 F 

. Es- 


timates of binary li fetimes, ^nsm 0.1-1 Gyr (e.g., 
iDominik et al.l l2012fl . are in conflict with the pres¬ 
ence of r-process-enhanced stars below [Fe/H] ^ —2.5 
in those models. Even if a short-lived channe l with 
/nsm 1 Myr exists (iBelczvnski fc Kalogeral 120011 : 
[Dc Ponder fc Vanbeverenll2004). the low Gal actic rate of 
such events (0.4-77.4 Myr~^. rPominik et al.ll20T^ leads 
to a too large star-to-star scatter in [E u/Fe] at [Fe/H] 
> —2.5 (lQia.nl I 2 OOOI : lArgast et al.l 1200411 . It should be 
noted, however, tha t the inhomogeneous GCE model of 
lArgast et al.l (120041 also appears to be in conflict with 
the observed small scatter of a elements relative to iron 
([Argast et al.ll200§ l. This problem would be cured if the 

^ N ote that the higher resolution model of Ivan de Voort et al.l 
lf20T^) also suggests short lifetimes of neutron star binaries. 
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TABLE 1 

Models of Sub-halos for Cases 1,2 


TABLE 2 

Parameters Adopted for NSMs 


M* (Mq) 

V 

koF (Gyr 1) 
Case 1 Case 2 

fcsF (Gyr 1) 
Case 1 Case 2 

10^ 

79 

1.0 

16 

0.013 

0.20 

10® 

40 

1.0 

7.9 

0.025 

0.20 

10® 

20 

1.0 

4.0 

0.050 

0.20 

107 

10 

1.0 

2.0 

0.10 

0.20 

10® 

5.0 

1.0 

1.0 

0.20 

0.20 

2 X 10® 

4.1 

1.0 

0.81 

0.25 

0.20 

inhomogeneity 

of the inter 

-stellar 

medium 

(ISM) 


modest because of effective matter mi xing (or if more 
mass ive CCSNe had greater iron yields. IWehmever et alJ 

It should be noted that the aforementio ned GCE mod¬ 
els as s ume a single halo system (except for lKomiva et al.l 
[Ml I van de Voort et ^ I2015I1 . which is incompatible 
with the currently dominant paradigm o f the hiera r chical 
merg ing scenario for galaxy formation. iPrantzo^ (I2006L 
I2008f) has shown that the overall shape of the Galactic 
halo’s metallicity distribution (MD) can be well repro¬ 
duced by a merging sub-halo model, with sma ller sub¬ 
halos having suffered larger outflows (see also iKomival 
201111. More importantly for the issue of r-process, 


PrantzosI (1200611 showed that if the sub-halos evolved at 


different rates, there would be no more a unique relation 
between time and metallicity; in that case, he argued 
that the observed “early” (in terms of metallicity) ap¬ 
pearance of r-elements and their large dispersion can be 
explained, even if the main source of those elements is 
NSMs. 

In this Letter, we study the role of NSMs in the early 
Galaxy with a GCE model based on the framework of 
such a hierarchical merging scenario and on recent esti¬ 
mates of the progenitor binary lifetimes. 

2. CHEMICAL EVOLUTION OF SUB-HALOS 

Each sub-halo is modeled as a one-zone, well-mixed 
single system that is losing gas by outflow (because 
of stellar winds, SN explosions and tidal interac¬ 
tions). The GCE model used in this work is the 
sa me as that of the ISM evolu ti on of the Galac t ic hal o 
in llshimaru fc Waiiaiol (1199911 : llshimaru et al.l (I2004D : 
IWanaio fc Ishimarul ( 2006[) but it differs in what follows. 
Each sub-halo is composed of stars with final total stel¬ 
lar mass M*, evolving homogeneously (i.e., with its ISM 
well mixed at every time). We consider sub-halos with 
M*/Mq = I0^ 10^ 10®, 10^ 10®, and 2 x 10®. The heavi¬ 
est mass is set to be about half the estimated stellar mass 
of the Galactic halo, ^ 4 x 1 O®M 0 (iBell et al.ll2008f) . The 
lowest mass corresponds t o the one of recent ly discovered 
ultrafaint dwarf galaxies (IKirbv et al.ll 200 ^ . 

Although we have no direct observational clues for the 
property of each sub-halo, recent observations of MDs 
of dw arf galaxies in the Local Group can provide some 
hints (iHelmi et al.l 1200^ . It is known that [Fe/H] at 
the peak of a MD, [Fe/H]peak, corresponds to effective 
yield yeff, i-e., the iron productivity of a galaxy system 
which depends on the adopted st ellar yields and IMF 
((Helmi et al.l 120061 : iPrantzod 1200811 . If a dwarf galaxy 
loses a significant amount of iron through outflow with 
the rate (OFR) proportional to the star formation rate 


Type 

iNSM (Myr) 

/nsm//ccsn 

Meu.NSM (Mq) 

short-lived 

1.00 

1.0 X lO-'^ 

2 X 10-® 

long-lived 

100 

1.9 X 10“® 

2 X 10"® 


(SFR), T/eff is ap proximately pro portional to r]~^, where 
OFR = 77 SFR (|Praiitzoi [2(30811 . In practice, limited 
data of MDs are available and thus the observed me¬ 
dian metallicities ([Fe/H]o) of dwarf galaxies are used 
instead of [Fe/H]peak- Observed mass-metallicity rela¬ 
tion of dwarf galaxies suggests that [Fe /Hjp approxi¬ 
mately scales as M®-® (IKirbv et al.ll201,'lll . We thus as¬ 
sume y~f^ oc 77 = 5.0 (M*/ 1 O®M 0 )“® ® as in iPrantzoi 
(Ml 0 . This can be naturally interpreted as a re¬ 
sult of smaller SFR and/or greater OFR because of the 
shallower gravitational well for a less massive sub-halo 
(see also §4). SFR is assumed to be proportional to the 
mass of ISM (Mjsm): SFR = /csfAIism, where the con¬ 
stant fcsF corresponds to the star formation efficiency, 
and OFR = kopMisM as well, where fcoF = vksp- How¬ 
ever, we cannot constrain either of fcsF or fcoF from ob¬ 
servations. Thus, we assume two extremes: Gases I and 
2 for the fixed ratios of i], respectively (Table 1). The 
sub-halo of M* = 10 ®M 0 is (arbitrary) set to have the 
same values of ksp and kop for both Gases 1 and 2. 

The GCSN abundances of Mg and Fe (which we as¬ 
sume to be representative of GGS N products and meta l- 
licity, respectively) are taken from lNomoto et all (l2006f) . 
We ignore the Fe contribution from Type la SNe, be¬ 
cause they appear to have played negligible role during 
the Milky Way halo evolution (the observed a/Fe ratio of 
halo stars is ^constant). The evolution of each sub-halo 
is computed up to 2 Gyr, which is considered a reason¬ 
able estimate of the time-scale of the Milky Way halo 
formatio n suggested from the observed ages of globular 
clusters (iRoediger et al.ll^014ll . The value of fcsF for the 
sub-halo with M, = 1 O®M 0 is set to account for the most 
metal rich stars in the Galactic halo. 

Binary population synthesis models indicate an av- 
erage binary lifetim e of (^nsm) 1 Gyr (e.g., 

iDominik et al.l l2012fl . being comparable with the life¬ 
times of sub-halos. A significant fraction of NSMs 
are, howe ver, expected to have tNSM 100 Myr 
(^ 40%, iBelczvnski et al.l 1200811 . In addition, a 
short-li ved channel of ^nsm 5, 1 Myr is also pre¬ 
dicte d (I Belczvnski fc Kalogeral l2?)f)ll: IBelczvnski et al. 
I 2 OO 4 [Pt^ Ponder fc VanbeverenI 120041 1. IDominik et al. 
l|2012ll estimate the fraction of such a putative short¬ 
lived channel to be up to ^ 7% (for the solar metal¬ 
licity case). For illustrative purposes, we assume here 
a bimodal distribution of ^nsm = 1 Myr and 100 Myr 
and with the corresponding fractions of 5% and 95%, 
respectively (Table 2). For the total average frequency 
of NSMs (/nsm) relative to that of GCSNe we assume 
(/nsm) = 2 X 10“®(/ccsn) adopting the e stimates of the 
recen t population synthesis calculations (|Dominik et al.l 
I 2 OI 2 I I. The mass of the ejected Eu, representative of 

^ The v alues of the coefficient and the index are updated from 
IPrantzol (120081 ') by adopting the latest mass-metallicity relation 
(IKirbv et al.ll2Qm . 
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TABLE 3 

Results Related to MDs of Sub-halos 


M* (Mq) 


Md {Mq) [Fe/H]peak 

Case 1 Case 2 Case 1 Case 2 



7.6 

X 

10® 

6.6 

X 

10® 

-2.63 

-2.56 

10® 

3.8 

X 

107 

3.3 

X 

10^ 

-2.33 

-2.30 

10® 

1.9 

X 

10® 

1.7 

X 

10® 

-2.03 

-2.02 

10^ 

1.0 

X 

10® 

9.0 

X 

10® 

-1.74 

-1.74 

10® 

5.3 

X 

10® 

5.3 

X 

10® 

-1.46 

-1.46 

2 X 10® 

8.8 

X 

10® 

9.3 

X 

10® 

-1.38 

-1.39 


r-process elements in the Galacti c halo, is taken from 
the latest nucleosynthesis result of lWanaio et al.l ( 20141) , 
Mf,ii nsm = 2 X 10~ ^Mr7^ (see also iGorielv et al. 120111 : 
IKorobkin et alJ[2012D . The mass of the ejected Ba, which 
is also predominantly produced by the r-process in the 
early Galaxy, is set to be -Mea.NSM = 9 Meu nsm ac- 
cording to the sol ar r-process ratio of these elements 
(|Burris et alJl200^ . 

3. ENRICHMENT OF r-PROCESS ELEMENTS 

Our results for Case 1 (with constant /cqf) are shown 
in Figure [T] Those results for less (more) massive sub¬ 
halos are indicated by the thinner (thicker) curves. The 
cumulative number of CCSNe in each sub-halo mono- 
tonically increases with time (Figure [T^). For NSMs, 
the cumulative number steeply rises when the long-lived 
binaries start contributing at 100 Myr. [Fe/H] mono- 
tonically increases with time for all the sub-halos; at a 
given time [Fe/H] is greater for more massive sub-halos 
because of their higher fcsF- MDs for all the sub-halos 
are also shown in Figure [TJ;. The metallicity at the peak 
of each MD, [Fe/H]peak (Table 3), is in agreement with 
the observed mass-metallicity relation (|Kirbv et al.l2013ll 
scaled downward by —0.4 dex to exclude SNe la contri¬ 
bution. The dark halo mass of each sub-halo, Md (Ta¬ 
ble 3), estimated from the initial baryonic mass and the 
initial baryon to dark mass ratio (which is assumed to be 
equal to the cosmic value Db/Dd = 0.046/0.24 ~ 0.19) 
implies Md oc Since the reasonable mass func- 

tion of Mp can be regarded as dTV/dMp oc M/(^ (e.g., 
iPrantzoi 120081 and the references therein), we obtain 
the sub-halo mass function as oc We find that 

the total MD (thick-black curve) weighted with the sub¬ 
halo mass function is in reasonable agreement wit h the 
observed one (gray-hatched region. lAn et al.ir2013D . We 
also find that the evolution of Mg (Figure[TJi, representa¬ 
tive of GGSN products) is in reasonable agreement with 
the observed stellar abundances of Galactic halo stars. 
In contrast to r-process elements, the observed scatters 
in [a/Fe] is known to be as s mall as the measurement 
errors le.g.. iCavrel et al.ll2004D . This result is consistent 
with such small scatters in [Mg/Fe], because each evolu¬ 
tionary trend is almost independent of M*. 

The resulting evolution of Fu (representative of r- 
process elements) is presented in Figure[T^ and compared 
to observed stellar values. We notice a transition from 
slow to rapid evolution of [Fu/Fe] for each sub-halo at 
~100 Myr, as a result of the contribution from long- 
lived binaries starting. The corresponding [Fe/H] (see 
Figure [TJj) differs from one sub-halo to another, being 
lower than [Fe/H] ^ —3 for M* < Mq. This in¬ 


dicates that the presence of stars at [Fe/H] ^ —3 with 
star-to-star scatter in [Eu/Fe] (< 0.5) can be interpreted 
as a result of NSM activity in sub-halos with various fcsF- 
Our model cannot, however, explain the presence of r- 
process-enhanced stars with [Eu/Fe] >1. In addition 
our result predicts the presence of stars having [Eu/Fe] 
^ —1 for [Fe/H] < —3. This is due to the contribu¬ 
tion of the short-lived NSM at early times (< 100 Myr). 
Measurements of Eu in such stars would be challenging 
because of the weak spectral lines. Such a signature has 
been seen in the Ba abundances of EMP stars, [Ba/Fe] 

^ —2-1 at [Fe/H] < —3 (Figure[p'), which could also 

be explained as due to the contribution of short-lived bi¬ 
naries. Note that the [Ba/Fe] values for [Fe/H]> —2.5 
are underpredicted compared to the observed ones; in 
this metalli city, contribution f rom the s-process becomes 
important ([Burris et al.ll2000|) . 

Figure [2] shows the results for Gase 2, where a con¬ 
stant /csF is assumed. We find the cumulative numbers 
of NSMs and GGSNe being similar to those in Case 1 
(Figure [T^). However, the evolutions stop earlier for less 
massive sub-halos. This is due to the termination of star 
formation because of gas removal by their significant out¬ 
flows with greater values of fcoF- In contrast to Case 1, 
the [Fe/H] evolution (Figure Od) is identical among sub¬ 
halos because of the same fcsF, although the termination 
points differ depending on their fcoF ■ The resulting MDs 
(Figure [2t) as well as the evolutions of Mg (Figure [2ji) 
are similar between Cases 1 and 2. In Figure [2fe, we find 
that the r-process abundances increase to values greater 
than [Eu/Fe] ^ 1 for the sub-halos with M* < 10® Mq. 
This is a consequence of the fact that the higher fcoF 
lead to smaller Fe amounts and thus higher Eu/Fe ra¬ 
tios. It is interesting to note that, even without invoking 
inhomogeneity of ISM, the enhancement of Eu can be in 
part explained by our sub-halo models. However, our re¬ 
sult here cannot account for the r-process enrichment at 
[Fe/H] ~ —3 because the [Eu/Fe]’s start rising at [Fe/H] 
~ — 2.4 for all the sub-halos. 

Our results imply that the reality may be between 
these two extremes, Cases 1 and 2; reasonable com¬ 
binations of ksF and /cqf might account for the presence 
of r-process-enhanced stars at [Fe/H] ^ —3. It is also 
important to note that the cumulative numbers for the 
least massive sub-halos (M* = 10"^ Mq) are ~ 0.1 around 
the end of their evolutions. This could be another source 
of large enhancements of [Eu/Fe] (> 1). The [Eu/Fe] 
values would be substantially higher than the averaged 
curves of our models, provided that only a fraction of 
sub-halos experienced NSMs. 

4. SUMMARY AND DISCUSSION 

We studied the role of NSMs for the chemical evolution 
of r-process elements in the framework of the Galactic 
halo formation from merging sub-halos. It is found that 
the appearance of Eu at [Fe/H] ~ —3 with star-to-star 
scatter in [Eu/Fe] (< 0.5) at [Fe/H] ^ —3 can be inter¬ 
preted as a result of lower star formation efficiency fcsF 
for less massive sub-halos. On the other hand, the pres¬ 
ence of highly r-process-enhanced EMP stars ([Eu/Fe] 
> 0.5) can be explained if values of fcoF (the multiplica¬ 
tive factor for the outflow rate) are higher for less mas¬ 
sive sub-halos. These may be reasonable assumptions be¬ 
cause less massive sub-halo systems have weaker gravita- 
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Fig. 1. — Evolution of the sub-halos with M^/Mq = 10^, 10^, 10®, lO"^, 10®, and 2 x 10®, respectively, indicated by the thinnest to thickest 
curves for Case 1. (a) Cumulative numbers of NSMs (solid) and CCSNe (dashed) as functions of time. The horizontal dashed line marks the 
number of unity (see text for implications), (b) [Fe/H] temporal evolutions, (c) MDs of sub-halos weighted with the sub-halo mass function 
and their sum ( thick-black). Observational data of the Galactic halo (gray-hatched histogram) are taken from the calibration catalog of 
lAn et akl J2Q13) . (d)-(f) [Mg/Fe], [Eu/Fe], and [Ba/Fe] as functions of [Fe/H], respec tively. The horizo ntal and vertical lines indicate the 
solar values. Observational stellar values (dots) are taken from the SAGA database IjSuda et al.ll20Q'^ . excluding carbon-enhanced stars 
that may have been affected by gas transfer in binaries. 







Fig. 2.— Same as Figure ^ but for Case 2. 

tional potential (as indicated by Md in Table 3) and thus 
are expected to form stars less efficiently and/or expel the 
ISM more easily. The ratio of OFR to SFR, 77 , is assumed 
to be proportional to Recent observations of rel¬ 

atively massive galaxies ( 10 ^ —IO^^Mq) also suggest sim¬ 
ilar a nti-correlation between M* and 77 (|Chisholm et al.l 
l2014fl . Under this assumption, the metallicity at the 
peak of each sub-halo’s MD (Table 3) appears to be con¬ 
sistent with the observed mass-metallicity relation. In 


addition, the total MD weighted with the sub-halo mass 
function shows reasonable agreement with the observed 
one of the Galactic halo. 

The low-level Ba abundances ([Ba/Fe] ^ —1.5) ob¬ 
served in the EMP stars of [Fe/H] < —3 may be due 
to the contribution from the short-lived binaries (with 
Gsm = 1 Myr in this study). Thus, the main features 
of r-process abundances observed in EMP stars, their 
appearance at [Fe/H] ~ —3 with large star-to-star scat- 
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ter and their sub-solar amounts for [Fe/H] < —3, can 
be potentially explained by the contribution from solely 
NSMs. 

The highly Eu enhanced stars at [Fe/H]~ —3 may be 
accounted by reasonable combinations of A:sf and fcoF- 
The small cumulative numbers of NSMs (< 1) for least- 
massive sub-halos could be additional reason for the ob¬ 
served large scatter in [Eu/Fe], because NSMs may oc¬ 
cur only in some of them, while no NSMs occur in others. 
Once NSMs occur in such less massive systems, their ISM 
must be highly enriched by r-process elements. 

In conclusion, our results imply that the current obser¬ 
vational evidences of r-process signatures in EMP stars 
can be interpreted as a consequence of the Galactic halo 
formation from merging sub-halos. If the model is a 
reasonable simplification for such chemo-dynamical evo¬ 
lutions of sub-halos, NSMs can be the main r-process 
contributors throughout the Galactic history. Contribu¬ 
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tions from CCSNe are not necessarily needed as invoked 
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